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Abstract 

In this short note we will use the residue theorem to establish a formula for Eu- 
ler's constant. In particular, we offer a slightly generalized version of an interesting 
infinite series due to Flajolet, Gourdon, and Dumas. 

1 Properties of 7 

Recall that 7 is usually defined in terms of the partial sums of the harmonic series: 



7 : 
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lim (V--logn I = 0.577215664901.... (1.1) 



While this definition is interesting, it is extremely inefficient for actually calculating 
the constant numerically. Since Euler's constant arises in many mathematical con- 
texts [5], and since it is often considered to be the most important "special constant" 
after ir and e, it is naturally an interesting problem to find rapidly converging series 
expansions for 7. 

While TT and e were both proved to be irrational in the eighteenth century, 
the arithmetic nature of 7 remains a mystery. Recall that Fourier obtained the 
irrationality of e as a simple consequence of Euler's famous infinite series: 

00 

k=0 

While it seems virtually certain that no identities like (II. 2p exist for 7, many inter- 
esting formulas are scattered throughout the literature. For example, Sondow and 
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Zudilin proved several rational expansions for 7 in [4j. A variety of series expansions 
involving the Riemann zeta function, and multiple identities containing logarithms 
also exist (see [2] or [5]). In the following theorem we offer a slightly generalized 
version of a formula due to Flajolet, Gourdon and Dumas (see page 28 in [1]). 

Theorem 1.1. Suppose that x > and w > {), then the following identity is true: 
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(1.3) 
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Proof. Let 4>w{x) be defined by 
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(t>w{x) ■■= y^e" 



k=0 

and suppose that w > 0. Since this sum converges uniformly it is easy to calculate 
the Mellin transform of (j)yj{x): 







Inverting the Mellin integral yields 



^Ux) = f^da, (1.4) 

where C is a closed contour which runs vertically along the line (^ — ioo, ^ + ioo), 
and then encircles the negative half plane. Notice that the integrand in (jl.4p has 
simple poles on the imaginary axis at a = 2Trik/w for k E {±1,±2, ...}, simple 
poles at the negative integers, and a double pole at a = 0. It follows easily from the 
residue theorem that 

„ , r(a)x~"\ ^ „ /T(a)x-''\ ^„ /T(a)x-'' 
(x) = Res 



a=o V 1 - e-"-"^ 



Res flWillVy-R.sfl 
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1 7 log(x) 1 ^ ^ /2vrifc\ , ~ (-1)'=+! x' 



+ 1 y rf^^x-2-^/-+y. 



w W W ^ \ w ) ' k\ e^^ — 1 
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Rearranging the last identity completes the proof. ■ 

Notice that a wide variety of formulas involving 7, vr, and values of the Riemann 
zeta function at odd integ ers can be obtciined by rnodifying the definition of (x) 
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(for instance by adding arithmetic functions to the sum). As a simple example we 
can show that 



CXI 



672 + ^2 ^ 1 + 12 ^ A: ( e""' - e~^'' + 1 ) - 24 ^ Re [r'(27rifc)] . 



k=l k=l 



We will also point out that equation (jl.3p reduces to a classical formula involving 
the exponential integral when w approaches zero [2]: 
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log(x)-y e-"d. + ^ ^ ^, 
•^•J k=i 

log(x) — / du + / du 
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To prove this limiting case, simply observe that the vanishing of the Gamma sum 
follows immediately from a standard estimate [3]: 



|r(27riA;/u;) | = O (^e^^ '^/'"^ 

as — 7- CX3. Finally, we can recover the formula of Flajolet, Gourdon and Dumas [l] 
by setting w = ln(2) in equation ()1.3p . 

Corollary 1.2. The following identities are true: 



fc=0 k=l ^ ' k=l 

fc=0 A:=l ■"" " ' k=l 

Proof. Equation ()1.5p follows from setting x = e and t/; = 2 in ()1.3p . while ()1.6p 
follows from setting x = w = l.B 

Despite converging rapidly (notice that Re [r(7riA;)] = 0{e~'^'^^^) as A; — )■ 00) 
both equations (jl.Sp and (II. 6p suffer from a serious computational drawback. In 
particular, there is no obvious way to calculate r(7riA;) recursively with respect to k. 
This fact most likely precludes the possibility of using either (jl.Sp or (II. 6p to obtain 
more than a few thousand decimal digits of 7. As of December 2006, the value of 
7 was known to more than 116 million digits [6]. Although a computational record 
may therefore seem far out of reach, 1 will point out that the first sum in equation 
(jl.Sp would only need around 9 terms to reach that level of accuracy (the tenth 
term is e'^^^ ^7 x 10-572,754,397 ^ ;l0-ii6'°°°'°°°). Thus, it may be a worthwhile 
endeavor to try to accelerate the convergence of the the second and third sums in 
both equations (II. 5p and (II. 6p . 
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